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Hybrid LES/subscale modeling approaches have an important advantage over 
the current noise prediction methods in that they only involve modeling of the 
relatively universal subscale motion and not the configuration dependent larger 
scale turbulence . Previous hybrid approaches use approximate statistical 
techniques or extrapolation methods to obtain the requisite information about the 
sub-filter scale motion. An alternative approach would be to adopt the modeling 
techniques used in the current noise prediction methods and determine the 
unknown stresses from experimental data. The present paper derives an 
equation for predicting the sub scale sound from information that can be obtained 
with currently available experimental procedures. The resulting prediction method 
would then be intermediate between the current noise prediction codes and 
previously proposed hybrid techniques. 


1. Introduction 

Crighton (1979), in a now classical paper, makes two major criticisms of the Liiley 
(1374) equation approach- which effectively treats the sound propagation as a 
small perturbation of a sheared mean flow. The first of these is that it provides no 
rational means of excluding the instability waves that eventually become infinite 
in any causal solution to the problem. This issue was eventually resolved by 
Goldstein and Leib (2005) and Goldstein and Handler (2003).The other criticism, 
which is dealt with in the present paper, is that the evolution time of the turbulent 
eddies is much greater than the time it takes a sound wave to traverse them- 
implying that the turbulence is effectively frozen during the passage of the sound 
and that there is, therefore, no particular mean flow profile that would give a true 
representation of the refraction effects that occur in the real flow. 

This latter difficulty would not occur if the sound field were determined from a 
DNS (Freund, 2001) or LES solution. Unfortunately, the computer resources 
needed to resolve all the relevant length scales are enormous — often well 
beyond those available on present day machines. A reasonable compromise 
would be to use LES simulations with very broad filter widths. This approach was 


adopted by Bastin, Lafon and Candel (1997), Bodony and Lele (2002, 2003, 
2005), and others. They found that simulations of this type are quite accurate at 
low frequencies but significantly under predict the high frequency component of 
the spectrum — presumably because they do not account for the sound generated 
by the sub-filter scale turbulence. Bodony and Lele (2002 ), attempt to resolve 
these difficulties by adopting a “hybrid” approach that models the non-linear 
Reynolds stress terms in an appropriate equation for the sub-scale motion and 
then uses the result to calculate the missing sub-scale sound — which can be 
added to the LES sound field. The great attraction of this approach is that it 
directly calculates the configuration dependent component of the sound while 
modeling only the relatively universal small scale motion. 

The relevant equations are obtained by dividing each of the flow variables in the 
Navier-Stokes equation into a component that is determined from the LES 
solution and a residual component that satisfies the Navier-Stokes equations with 
the LES equations subtracted out. By introducing new (in general non-linear) 
dependent variables (Goldstein, 2002, 2003) the latter equations can be rewritten 
in the form of the linearized Euler equations with sources that are the same as 
those that would be produced by external stress and heat flux perturbations. The 
corresponding source strengths, which depend on the non-linear sub-scale 
stresses, can, in principle, be modeled and the resulting linear equations can 
then be solved by using a Greens’ function approach. The result can then be 
used to obtain an expression for the far field pressure autocovariance in terms of 
the sub-scale turbulent motion. Unfortunately it is very difficult to model the 
instantaneous subscale stresses that appear in this result. 

Bodony (2004), in his outstanding Ph.D. thesis, uses an approximate 
deconvolution method to kinematically extrapolate information about the filtered 
scales to obtain the required information about the sub-filter scale motion. But, as 
Bodony (2004) points out, this is still too expensive to be “practical in an 
industrial noise prediction setting”. A more practical approach might be to extend 
the modeling approach used in current noise prediction methods such as the 
MGB (Balsa, Gliebe, Kantola, Mani, Strings and Wong, 1978) and JeNo 
(Khavaran, and Bridges, 2004) codes that determine the unknown stresses from 
experimental data. The hope is that the sub-filter stresses will be much more 
universal than the large scale stresses that have to be modeled with these 
current methods and that models developed from any flow will apply to any other. 
This approach would be intermediate between the existing hybrid approaches of 
Bastin, et al.(1997), Bodony and Lele (2003) and Bodony (2004) and the current 
(i.e., MGB and JeNo) noise prediction methods. 

But most experimental reports only document the reproducible (i.e., non- 
random) characteristics of these sources, such as their lower order statistics. It is 
therefore important to derive a formula for the subscale sound field that depends 
only on non-random quantities the can be measured with currently available 
experimental techniques. Bodony (2004) and Bodony and Lele (2003) obtain 
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such an equation by dividing the flow into steady and fluctuating components, 
deriving an inhomogeneous wave equation for the latter and then further sub- 
dividing the unsteady Reynolds stress source terms in that equation into filtered 
and subscale components- with the latter part being inputted through an 
appropriate model. Unfortunately this approach does not explicitly account for the 
scattering of the subscale sound by the filtered scale turbulence and therefore 
does not overcome Crighton’s (1979) second criticism of the Lilley equation 
approach. This is analogous to the arguments that eventually led to the 
replacement of the original Lighthill (1952,1954) approach with the Lilley (1974) 
equation formalism. While Lighthill's equation is exact with all real flow effects, 
including the mean flow interaction effects, contained in the source term, it tuned 
out to be almost impossible to model these weak but non-local effects from 
experimental measurements. This should also true for the non-local scattering 
effects and the same arguments would dictate that they be included as part of 
the propagation effects. 

The present paper (which is base on the conceptually simpler double 
decomposition formulation) is an attempt to eliminate this problem by exploiting 
the statistically independence (i.e., the decoupling) of the filtered motion in an 
LES simulation from the detailed sub-scale fluctuations. The only coupling that 
can occur in these simulations is through the filtered Reynolds stress — which is 
eventually modeled in terms of the filtered variables and their derivatives. The 
LES solution is therefore calculated from a closed set of equations that only 
involve the filtered variables. In the “hybrid” approach, the deviation of the sub- 
scale stresses from the filtered stresses, say cjU , also has to be modeled. The 

present approach assumes that the model will be constructed from spatially 
filtered experimental data (which can be obtained from PIV measurements). 
Unlike the LES simulation, the instantaneous random subscale fluctuations in the 
experiment will eventually produce 0(1 ) changes in the large scale motion (due 
to the chaotic nature of the solutions to the N-S equations) and the entire 
experimental flow should therefore be uncorrelated (in the statistical sense that 
covariances of various quantities will vanish) with the LES simulation where the 
randomness comes from the random initial conditions-- which are themselves 
uncorrelated with the initial conditions in the experiment (Lesier, 1987, p. 230 and 
Pope, 2000, pp. 612-613). It therefore seems reasonable to model crL as if it were 

statistically independent of the LES solution. Notice that this does not imply that 
the expectation or average value of this quantity should be independent of the 
LES solution. 

The main purpose of this note is to show that this type of modeling leads to an 
expression for the far field pressure autocovariance that depends only on the 
overall correlation of the subscale stresses, which is a non-random quantity of 
the type that is usually reported in the experimental literature and not on their 
instantaneous values, which is not. The present approach also provides a 
framework for the systematic introduction of additional approximations to further 
simplify the results — with the most significant of these arising from the near 
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isotropy of the subscale turbulence. The most general subscale turbulence 
correlation tensor has 45 independent components, but the isotropy assumption 
reduces that to 7. And if retarded time variations (i.e., variations of the 
propagation factor over the correlation volume of the turbulence) are also 
neglected, only two of those actually enter the formula for the far field pressure 
autocovariance. it would be very difficult to argue for the neglect of retarded time 
effects if the non local effects had to be included as part of the source term. 


The relevant sub-filter scale acoustic equations are introduced in section 2 and a 
formal Greens’ function solution is obtained in section 3. This solution is then 
used (in section 4) to derive a formula for the acoustic pressure autocovariance 
and the statistical independence assumption/model is then used to eliminate the 
dependence on the instantaneous subscale fluctuations and thereby obtain a 
formula that involves only commonly reported statistical quantities. The result is 
simplified by neglecting retarded time variations and an isotropic turbulence 
model is introduced in section 5 to simplify it even further. 


The pressure p , density p and velocity Vj in any flow can be decomposed into 

filtered p ,p and Favre (1969) filtered v- /p components that satisfy the LES 
equations and residual components 

p' = p-p, p'^p-p C 1 - 1 ) 

that satisfy the 5 formally linear equations (Goldstein, 2002, 2003) 


Z^ v w v = s p for \x, v = 1, 2, 3, 4, 5 
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and 

h 0 =h + -V 2 , %mh’ + -v' 2 (1.11) 

2 2 

where c 2 -y Rf is the Favre filtered square sound speed and, as usual, the 
viscous terms, which are believed to play an insignificant role in the sound 
generation process, have been neglected. 

3. Vector Greens 5 function solution 

These equations can be solved in terms of the vector Greens’ function 
g va (x,r|j,x) which satisfies the inhomogeneous linear equations 

(2- 1 ) 
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with delta-function type source term together with the causality condition 

g = 0 for t< x 


( 2 . 2 ) 
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From the acoustics perspective the primary interest is in the 4 th (Le., the 
pressure-like) component of (3.3) which only involves the 4 th component vector 

Greens’ function^ {x,t\y,z) . The latter quantity can, in principle, be found by 

solving be solving the system (3.1). But since this consists of 15 first order 
equations, it turns out to be simpler to compute the adjoint Greens’ function 

g a v 4 (y,T\x,t) , which was introduced into the Aeroacoustscs literature by Tam and 

Auriault (1998) to study sound propagation through steady flows and used 
effectively by Bodony (2004) to implement his hybrid approach. Equations (A-1) 

to (A-4) show that it is related to g 4v (x,t\y t x) by the reciprocity relation 


g a v4 ( y^W) = s 4 v { x ^\y> T ) 


(2.5) 


and satisfies the system 


) jr g ° 4 {y^ I-m) = - y)s(t-r ) 


( 2 . 6 ) 


of five first order equations which are given in component form by (A-5a,b,c). 
These equations suggest that g v ^{y\x:x,t-x) = g^^y^x.t^ will be a stationary 

random function of its third argument, x , which results from the time dependence 
of the coefficients in that are determined from the LES solution and are, 
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therefore, random functions of time. (Notice that this argument would not be 
needed if the LES solution were steady.) 

It therefore follows from (2.7), (2.8), (3.4) and (3.5) that the 4th component of 
(3.3) can be written as 

oo 

= J Y y [y\x\x,t-x}ay{y,x)dydx (2.7) 

V -00 


where like, g v n(y\x:T.,t -t) = gl t {y,x\x,t) , the propagation factor 
along with 
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of X jk (j|x:T,r-x) should be random functions of their third arguments (which 
arise from the third argument of g vA {y\x:x,t-x) as well as the x -dependence of 

dv k /dyj). 



Since n '(x,t) — -> p'[x,t) asx — » oo , the far field pressure autocovariance 




lim 

x — >00 


lim — 

r^oo 27 


i 

| p' (xj)p' (x,t + x)dt 


(3-1) 


becomes 


n=-i- } | x:t’,t’,T + t’-t’) +>/, t’)dy dtjdt' dt" 

LI 


T oo 


_i_ 

2r 


-T-o oV 


T oo 


(3.2) 




(y,fj | jt:r',?"+/',T + r") c^- (j/) (j + ^,t" + t')dydtjdt' dt n 


-T- oo F 


where 


(3.3) 

oo 


is a stationary random function of its fourth and fifth arguments. 

T iaJil (yJi | *:r',fW,T + ;")and 

R o n y(y’ , }’ t "’ t ')-^i O’.'') +t l’ t "+ t ' ) (3-4) 

are, therefore, both stationary random functions of t' . They should be considered 
to be uncorrelated in the present approach because, as note in the Introduction, 
the filtered motion, which is determined from an LES solution to a closed set of 
equations that only involve the filtered variables, should be statistically 
independent of (i.e., decoupled from) the detailed sub-scale motion determined 
from actual experimental data because, unlike the LES simulation, the 
experimental subscale fluctuations eventually produce 0(1) changes in the large 
scale motion (due to the chaotic nature of the solutions to the N-S equations). 

This causes the entire experimental flow to be uncorrelated with the LES 
simulation where the randomness comes from the random initial conditions-- 
which are themselves uncorrelated with the initial conditions in the experiment 
(Lesier, 1987, p. 230 and Pope, 2000, pp.812-613). The only coupling that can 



occur is through the filtered Reynolds stress — which is eventually modeled in the 
LES simulation. It is therefore reasonable to require that any model used to 

represent the deviation cr' u/ =-(pv^vt + <? H/ j of the sub-filter Reynolds stresses 
from the filtered stresses & U7 be statistically independent of the filtered motion 
(and consequently with the coefficients ofZ,”„ ). This means that the covariance 
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over the product of two non-random functions. 
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r ,*(.F,»/ | x\t",z+t") = J y m (y# | x:t",t + z + t’,t) dt (3.10) 


over time of the correlation 

1 T 

y m {yt! I x:t',t+x + f,t) -//m — |y te (j'|-*::r',/ + T + <")y, M + + (3.11) 

of the two random propagation factors y,„ (y\x.t',t + x + t"-t') and y ((l (j + jj|x:r',r) , 

can be thought of as an expected or mean propagation factor. It can be 
determined as part of the LES computation. The second factor, which, can be 
thought of as a source function, describes the low order statistics of the unknown 
component of the sub-scale fluctuations. It has to be modeled, but it should be 
much easier to do this than to model the instantaneous values of the sub-scale 
fluctuations (which would have to be done if the source and propagation 
fluctuations were not de-correlated). Notice that this de-correlation only implies 
that the fluctuation, F" a/> , (and not T m itself) be independent of the resolved 

scales. It may, therefore, be appropriate to paramatize the model for this quantity 
and determine the parameters from the LES calculation. 

It is useful to introduce the pressure autocovariance 


00 

u(x\ y,x) = J \f ta] ly,n I x :?,* + ?) R oi ^{y ;n ,t") di]dt" asx^ CO (3.12) 
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at the observation point x due to a unit volume of subscale turbulence at the 
source point y- which only involves the statistical correlation R G¥J {y ,t n ) of the 
subscale velocities and not their instantaneous values. 

But computation of the far field pressure still requires a great deal of information 
about these quantities (which cannot be obtained from the LES solution). This 
requirement would certainly be reduced if the correlation volume of the subscale 
turbulence were small compared to the characteristic length scale over which the 
filtered scale turbulence, and therefore the propagation factor 
TiJy?l | x:t",T+t' r ^ , varied. The latter could then be treated as a constant 
relative to the n -integration causing the result to depend only on the integral of 
R cr ^(y ,f) over the separation variable x \ , which, in turn, depends only on the 
decay time f at any given source point y . Since most of the non-local 
scattering effects are already accounted for by the propagation factor, this is 
likely to be the case for variations perpendicular to the mean flow direction, but 
not necessarily in the in the mean flow direction itself. However, as Lighthill 
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(1952, 1954) pointed out in a slightly different context, the streamwise decay rate 
should be much more rapid in a reference frame moving with the subscale 
convection velocity, UJ y ) , of the turbulence. Ffowcs Williams (1963) showed 
that this idea can best be implemented by introducing a moving frame 
correlation, say 




(3.13) 


where 


% = n~U c {y)x (3.14) 

before integrating the turbulence correlation over the separation vector. 
Introducing this into (4.12) and changing the integration variable to £ and 

neglecting variation of I] . over the correlation volume yields 


n(x|^,x)= \T m (y,Uf 
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Enthalpy fluctuations are thought to be negligible in cold air-jets (Lslley, 1 996). We 
therefore neglect these quantities, for simplicity. In which case (3.15) becomes 


IT(.x|jm)= \T m (y,U/\ + 

— oo 

as x — » oo 

Bodony’s (2004) LES computations imply that the convection velocity can be 
expected to vary with the source location y and to be scale dependent for the 
larger turbulent scales, but should become much more constant and scale 
independent for the smaller scales. We therefore expect the subscale convection 
velocity UJy ) to exhibit similar behavior. 

Obviously 
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(3.17) 


n(x,x)= \tl(x\y,x)dy 
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The 81 component tensor RJ ilk {y ; £ ,/") has 45 independent components. 

Fortunately, the small scale turbulence should be quite isotropic in the moving 
frame (but, unlike the medium scales, not even approximately quasi-normal). 
Batchelor (1953) points out that the most general 4 th order isotropic tensor is of 
the form 
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where 
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The latter quantity and, therefore, the G k ,k = 1,2, can be computed from the six 
independent components of X jk ~~t) which are given in terms of the 4 th 

component adjoint vector Greens’ function g a v ^{y,r\x,t) by (3.10). 


A previously derived equation for the sound generation by the sub-filter 
turbulence scales in a large eddy simulation (Goldstein, 2002, 2003) is used to 
obtain a formula for the acoustic pressure autocovariance in terms of the sub- 
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scale turbulence correlation tensor by exploiting the statistical independence of 
the filtered and sub-filter scale motion. The sub-scale acoustic radiation can 
therefore be calculated by introducing appropriate models for this relatively 
universal correlation while determining the configuration dependent large scale 
sound directly from the LES solution. This approach is in accord with recent 
speculation that there are two distinct sources of jet noise — one of which is 
associated with the large scale motion and the other with the small scale motion 
(Tam, 1998). The approach also accounts for the scattering of the sub-scale 
sound by the large scale motion-an effect that was emphasized by Crighton 
(1979). 

Finally, I would like to thank Dr. Daniel Bodony for commenting on the 
manuscript and Dr. Jayanta Panda for showing how the subfilter correlations can 
be determined from PIV measurements. 


The adjoint Greens’ function g a va (y,?\x,t) satisfies the adjoint equation (Morse 
and Feshbach, 1953, Tam and Auriault, 1998, Bodony and Lele, 2002 ) 
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and is related to the direct Greens’ function by the reciprocity 


relation 
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When written out in full the system (3.6) becomes (Bodony and Lele, 2003 
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which shows that satisfies the inhomogeneous wave equation 
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The relevant solution which satisfies the causality condition 



IS 



1 


a 


4 rt x-y c~ dr 


T~t + 


\x-y\ 


'0 J 


So that when x = \x ->oo 


a 


&44 


4 7TXCq dr 


f X'V x ^ 

T 1 H 

xc 0 C 0 J 


(A- 10) 


Notice that the far field adjoint Greens’ function can be calculated directly by 
allowing x = |x| ->oo in (A-7). . Smoother boundary conditions can be obtained by 

putting 

g44{}’’ T \ x ^) = 4frxc 2 dr 2 G * ( y ’ T \ x,t ) (A-11) 

which satisfies 
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together with the smoother boundary condition 
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where H denotes the Heaviside unit function. It may also be convenient to divide 
the solution into incident and scattered components. 
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